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We use the low frequency (110-180 MHz) capabilities of the Westerbork Synthesis Radio Telescope (WSRT) to characterise a large 
collection of single pulses from three low magnetic field pulsars. Using the Pulsar Machine II (PuMa-II) to acquire and coherently 
dedisperse the pulsar signals, we examine whether the bright pulses observed in these pulsars are related to the classical giant pulse 
emission. Giant pulses are reported from PSR Bll 12+50 and bright pulses from the PSRs Bl 133+16 and B0031-07. These pulsars 
also exhibit large intensity modulations observed as rapid changes in the single pulse intensity. Evidence of global magnetospheric 
effects is provided by our detection of bright double pulses in PSRs B0031-07 and Bl 133+16. Using the multi-frequency observa- 
tions, we accurately determine the dispersion measures (4.844±0.002 for B 1 133+ 16 and 9.1750±0.0001 for B1112+50), derive the 
radio emission height in PSR B 1 133+16 and report on the properties of subpulse drift modes in these pulsars. We also find that these 
pulsars show a much larger intensity modulation at low sky frequencies resulting in narrow and bright emissions. 

Key words, pulsars - neutron stars - emission mechanism - giant pulses 
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cd 1. Introduction 



__ , Most radio pulsars are characterised by a stable average pulse 

■ profile at a given sky frequency. The average profiles are formed 
f — | by integrating all pulses emitted by the star over a given time, 
r*"** ■ The pulsed emission from the star is observed as individual 
OO ■ pulses that carry detailed information about the physics of the 

| pulsar radio emission. These pulses can only be directly detected 

■ in a fewer than a third of the known radio pulsars, mainly be- 
7— I ' cause of their weak nature. One of the remarkable features of the 
CD single pulse emission is the occurrence of intense pulses called 
^ — 1 ■ giant pulses, which are defined as pulses with energy greater 

| than 10 times the average pulse energy. Furthermore, the gi- 
• rH , ant pulses are very narrow compared to the width of average 

■ pulsed emission and result in a pulse intensity distribution de- 
H | scribed by a power law. For example, in the radio emission of 

. . . i the young Crab pulsar, the giant-pulse widths are found to be 
narrower than 0.4 ns, implying a high brightness temperature of 
~10 41 K (lHankins & Eilekjl2007l) . and the singl e pulse flux dis- 
tribut ion follows a power law with slope -3.3 dLundgren et al.l 
1 19951) . Another example is the emission of single bright pulses 
from the young Crab-like pulsar, PSR B0540-69 located in the 
Large Magellanic Cloud; the pulsar is visible at radio wave- 
lengths only by virtue of th e star's intense giant pulse emission 
(I Johnston & Romanll 120031) . Similar bright and narrow pulses 
were reported in the emissi on of very old mil lisecond pulsars 
such as PSR B1937+21 ([Cognard et alJ Il996h and B 182 1-24 
(iRomani & Joh nston 2001). Concerted giant pulse searches in 
these older and short period pulsars have revealed more gi- 
ant pulse sources e.g., PSRs J1823-3021A and J0218+4232 
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dKnight et alJl2005l [2006bh . The only common feature of these 
pulsars is a strong magnetic field (Bic ~10 5 G) at the velocity of 
light cylinder, which is the limiting radius where the velocity of 
the co-rotating plasma and magnetic field lines is less than the 
speed of light. 

However, some studies propose giant pulse-like emission 
in pulsars with a low magnetic field at the light cylinder e.g., 
Bic in the 4-100 G range. For instance, the bright single pulse 
emission from PSR s B0031-07, J1752+2359, and Bll 12+50 
at 40 and 11 1 MHz (iKuzmin & Ershovl2004tlErshov & Kuzminl 
2005L 120031) . iKramer et alJ d2003l) also report bright pulses from 
PSR B 1133+16 at 5 GHz. If these bright pulses are simi- 
lar to the classical giant pulses, this questions the requirement 
of a high value of Bic to produce giant radio pulse emis- 
sion. After the initial discovery, a detailed statistical analysis 
of the bright pulses in these low Bic pulsars has not been at- 
tempted, unlike the extensive studies of the gi ant pulses from 
the y oung pulsars and the millisecond pulsars {L undgren et al 



1 995 t|Kinkhabwala & Thorsettl2000t lKnightl2007LlKnight et al 
I2006al) . 

The pulsars in this study also exhibit other interesting 
single pulse behaviour: PSRs B1 133+16 B0031-07, and 
Bll 1 2 +50 show drifting s ubpulses dBacket 19731: iTavlor et al] 
1 1975 1 : IWrightet alJ 119861) and nulling (iHerfindal & Rankinl 
2007; Ritchingsl 1 1 9761: IWright et all 1 1 98rtTvivekanandl 1 1 9951) 
PSR B1133+1 6 also shows narrow emission fe a tures called 
micro structure dHankinslll972l : ICordes et alJll990t iLange et all 
119981) . 

Radio pulsar s in general show a negative spectrum 
dSlee et alJ Il986l) and a low frequency spectral turnover 
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dKuzmin et alJ Il978t iMalofeev etall 11994 . This implies that 
pulsars are brighter at low sky frequencies, although the sky tem- 
perature also increases, contributing significantly to the system 
temperature. Moreover, microstructure and drifting subpulses 
can be more easily studied at lo w frequencies whe n sufficient 
telescope sensitivity is available dCordes et alJI 199Qh . Motivated 
by the possible presence of giant pulses and a rich single pulse 
behaviour in the pulsars considered here and the availability of 
the Low Frequency Front Ends (LFFEs) at the WSRT, we un- 
dertook this study in the 110-180 MHz frequency range. The 
flexible baseband recorder PuMa-II (Karuppusa my et al.| [2008) 
and full coherent dedispersion, permitted a detailed single pulse 
study with high sensitivity and relatively high time resolution. 
Although these pulsars have been studied at low sky frequencies 
before, we are able with our system to observe for much longer 
and across a far wider frequency range. Our long observations 
are useful in determining the flux distribution and the occurrence 
rates of bright pulses. Furthermore, the wide frequency coverage 
permits the computation of the pulsar spectra in this band. 

The rest of the paper is organised as follows. The observa- 
tions, the pulsar fluxes, and spectra are described in Sects. [2]and 
[3] We present our revised values of the dispersion measure (DM) 
of two pulsars based on the narrow giant pulses in Sects. [5]and[4] 
The single pulse intensity, energy distributions, and microstruc- 
ture are the subjects of discussion in Sects. [6]and[7j We proceed 
with the pulse drift analysis in Sect. [8] which is followed by our 
discussion and conclusions. 

2. Observations and data reduction 

The pulsars were observed with the WSRT on three different 
days between June-November 2008 in the 1 10-180 MHz range 
(see Table Q]). The design of the WSRT backend systems al- 
lows the selection of a maximum of eight different radio bands 
within the frequency range of the front end receiver. Eight 2.5 
MHz-wide bands tuned to various sky frequencies (see Table [TJ 
were chosen for this study since these bands showed compara- 
tively little radio frequency interference (RFI). The signals were 
recorded and processed using the PuMa-II instrument. When 
possible, the observations were carried out in the early hours of 
the day to reduce the effect of RFI. Despite our careful choice 
of the frequency bands and observing times, these observations 
were still susceptible to RFI because of (1) strong in-band RFI in 
the passband of the LFFE electronic amplifiers (2) the presence 
of the automatic gain control (AGC) in the WSRT's signal chain 
and (3) the 2-bit system design. Therefore, the data were sub- 
jected to an extensive RFI cleaning procedure by a combination 
of automated algorithms and visual examination. 

The WSRT was operated in the tied-array mode by adding 
signals from 13 telescopes each 25-m in diameter. The analogue 
signals from the individual telescopes were digitally sampled at 
2-bit resolution and a fixed rate of 40 MHz. The 2-bit data were 
then coherently added in a dedicated hardware unit to produce 
the 6-bit tied-array data. Coherent addition improved the signal- 
to-noise ratio (S/N) by a factor of 13. The summed data was then 
read into PuMa-II as 8-bit numbers, downsampled by a factor of 
8 in realtime, and written to disk. Single pulses were generated 
from the coherently dedispersed baseband data using the disper- 
sion measure in the catalog an d based on a polynomial deter- 
mined by the TEMPO software dTavlor & Weisberglll989h . 

The open-source DSPSBQ software package was used to 
form a 64-channel software filterbank simultaneously at the 

1 http://dspsr.sourceforge.net/ 



dedispersion stage in each 2.5 MHz band, and the resulting 
frequency-time cubes were written to the disk as single pulses. 
The original time resolution of 2.5 fis and the formation of a 64- 
channel synthetic filterbank permitted 2048 samples to be ex- 
tracted across the pulse period of the three pulsars observed and 
resulted in the final time resolutions in Col. 9 of Table Q] The 
large number of frequency channels was later used to correct for 
the residual dispersion smearing, since the DMs in the catalog 
had either changed or were not sufficiently accurate and it also 
aided the removal of narrow band RFI in the signa l. Extensive 
use of the PSRCHIVE package dHotan et alj|2004l) was made 
in the computation of the total intensity, signal-to-noise ratio, 
and visual inspection of the single pulses for RFI. Pulse stacks 
were formed by stacking individual pulse intensities one above 
the other according to their phases. When required, the signal 
was flux calibrated using the off-pulse radiometer noise. With 
the computed fluxes in each band, various statistical analyses of 
the data was performed. 

3. Low frequency flux 

3. 1 . Radio frequency interference 

The data was cleaned of interference using various cleaning 
strategies for the three purposes of estimation of mean flux den- 
sity, fluctuation analysis, and single pulse flux estimates. For the 
mean pulse flux computation, the pulses that showed excessive 
noise, negative drop outs, and large intensity peaks in the off- 
pulse region were removed. The single pulses were combined 
to create pulse stacks that were then examined for the second 
time as time-phase plots; impulsive interference missed by the 
automated algorithm was then identified and removed. The re- 
maining pulses were then averaged over the period of integra- 
tion to compute the average pulse profile. The number of pulses 
removed in this method was accounted for in the time term of 
Eq. [T]before the profile was flux calibrated. From the calibrated 
profiles, the spectral indices were derived. 

For the fluctuation analysis, the pulses were first gated to 
contain at least three times the width of the average pulse pro- 
file at the lowest observed frequency. The gated pulse stack 
was then subjected to the two-pass RFI cleaning procedure ex- 
plained above. Pulse stacks were inspected by eye and pulses 
that showed a saturation in the off-pulse phases were replaced 
by zeros. This procedure was repeated until a visually clean im- 
age was seen in the time-phase plane. 

The single pulse statistics were computed using only the 
single pulses from one of the eight bands observed. RFI was 
cleaned from the dedispersed single pulses by visually exam- 
ining the frequency-phase images within 30° of the pulse lon- 
gitude. For these comparatively long period pulsars, retaining 
a narrow phase range of the pulse period has the advantage of 
being less sensitive to impulsive RFI affecting other pulse lon- 
gitudes. Only pulses that showed the required dispersion were 
retained. An interactive plot utility was developed to examine 
the frequency-phase images of the large S/N pulses. 

3.2. Flux calibration 

Calibrating pulsar fluxes at low sky frequencies at the WSRT is 
difficult for three reasons. First the main contribution to the sys- 
tem noise T sys , is from the sky temperature 7^,. which in turn 
depends on the position of the pulsar on the sky. Second, the 
aperture efficiency of the LFFE's at the WSRT is poorly con- 
strained and thirdly, the presence of RFI introduces errors into 
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Table 1. Observation details and telescope parameters. 



Pulsar 


Period 


Date 


Duration 


N pulses 


Band" 


T 

A sys 




Af b 


V c 

u mm 




s 


dd/mm/yy 


min 


MHz 


K 






Jy 


B1133+16 


1.187913065936 


23/11/08 


280 


-13701 


110-180 


400 


2048 


809 


9.0 


B0031-07 


0.942950994559 


02/11/08 


360 


-22758 


110-180 


400 


2048 


580 


10.7 


B1112+50 


1.656439759937 


01/10/08 


570 


-14286 


110-180 


400 


2048 


460 


12.0 



Notes. 

(a) The pulsar signal was recorded from eight 2.5 MHz bands centred at 116.75, 130, 139.75, 142.25, 147.50, 156, 163.50, and 173.75 MHz. 

(6) Af is the effective time resolution of the pulses with N w „, across the pulse period. 

(l) 5„„„ is the minimum detectable signal per time sample based on a 30% aperture efficiency. 



the flux measurements of the calibrators and the pulsars. The 
system parameters (telescope gain and system temperature) were 
derived from continuum observations of the calibrator source 
3C196 made shortly before these measurements. From these 
measurements, the aperture efficiency is estimated to be ~30% 
translating to an effective telescope collecting area of 1500 m 2 . A 
constant value T sys in the entire 1 10-180 MHz is used, since the 
Galactic synchrotron foreground varying with frequency as v~ 2 6 
is almost removed by the frequency-dependent amplitude gains 
of the LFFEs. The system temperatures derived in this manner 
are shown in Table. Q] Assuming a frequency scaling of v~ 2 - 6 , the 
sk y temperature estimated from the 408 MHz continuum maps 
of Haslam et al] (1 198 lh in the direction of these pulsars is con- 
sistent to within 20% of the system temperatures listed in Table 
ffl 

Once the telescope gain and system temperature had been 
established, the off-pulse noise in flux units was computed using 
the radiometer equation (Dicke 1946) that can be expressed as 



where S sys = T sys /G is the system equivalent flux density in Jy, 
N p the number of polarisations summed, and T- m , the equivalent 
integration time. The pulse profile was then flux calibrated by 
multiplying by S av . Integrating the fluxes in all phase bins and 
dividing by the total number of bins then resulted in the mean 
continuum flux density of the pulsar. While this method is not 
as accurate as the absolute flux measurements based on noise 
generators, this is the only option in a synthesis telescope such 
as the WSRT for non-imaging observations; it also yields fluxes 
comparable to those reported elsewhere. 

In our method of detecting the pulsars, the dedispersed signal 
was written to disk as single pulses. Following the formation of 
total intensities, diagnostic plots of peak S /N against the phase 
of the peaks were made. Any pulsed emission could then be seen 
as clustering of data points near the nominal emission phase and 
the corresponding pulses were visually examined. Combining 
the single pulses in time after removing the RFI afflicted pulses 
produced the average emission profiles. Thus, this method is sen- 
sitive to any pulsed emission above the noise floor. 

3.3. Profile evolution 

The average pulse profiles in all recorded bands for the three 
detected pulsars are displayed in Fig. Q] PSR Bl 112+50 dis- 
played 10% peak widths Wio = 35.5-38.8 ms in this fre- 
quenc y range. This can be compared to Wio = 35.0 ms at 405 
MHz (Lori mer et al. I [1991 . The single component in the aver- 
age emission of PSR Bl 1 12+50 at these frequencies and at 328 



MHz dWeltevrede et alj|2006al) evol ves to a two-compo nent av- 
erage emission profile at 1400 MHz dWright et al.lll986l) . 

PSR B00 3 1-07 shows W w = 120.2-155.1 ms, while 
lLorimer et"all d 19951) report a value of 104.5 ms at 405 MHz. 
The trailing edge in the average profile of PSR B003 1 -07 shows 
a gradual decrease in steepness with frequency. We rule out scat- 
tering in the ISM as a possible cause of this decrease because 
that would manifest itself as pulse broadening by a factor of 
-5 in the 110-180 MHz band considered here, which is clearly 
not seen (e.g Figure Q]). However, this change in slope is in- 
dicative of the emergence of a second emissio n component at 
even lower frequencies e.g at 40 and 62 MHz (llzvekova et al.l 
1993) and represents almost the opposite behaviour of the aver- 
age emission profile in PSR Bl 1 12+50. Another contribution to 
the slope of the trailing part of the average emission profile in 
PSR B003 1 -07 could be the vertically drifting subpulses, which 
is discussed in Sect. [8] The characteristic age of PSRs Bl 1 12+50 
and B0031-07 is 10 Myr and 36 Myr, respectively. The general 
expectation of a narrower pulse profile for older pulsars is not 
valid with respect to these pulsars, as the impact angle of the ob- 
server's line of sight to the emission beams from the pulsars may 
not be equal. 

PSR B1133+16 displays a profile wid th ofW m = 44 0-48. 1 
ms compared to 41.8 ms at 405 MHz (Lori mer et alj [l995h . 
The component separation of the average profile decreases with 
frequency and this phenomenon was successfully explained by 
the radius-to- frequency mapping in the pulsar magnetosphere 
(ICordesll 978). The average profile also shows a double emission 
component, with a pronounced bridge emission that rises from 
14% to 25% of the peak flux density at our lowest and highest 
frequencies, respectively. The second component in the average 
emission profile also reduces from 90% to 70% of the leading 
component intensity with increasing frequency. Th ese variations 
clearly probe different parts of the emission beam. iNowakowskl 
(1996) assumed a model in which the pulsar beam consists of 
two concentric emission cones to explain the observed two- 
component emission profile and the increasing bridge emission 
with frequency. However, this model is unable to explain why 
multiple emission components are not visible at higher frequen- 
cies. Polarisation observations of this pulsar over a range of 
frequencies would offer further insights into th e emission ge- 
ometry . As discussed later in this work and by Kard ashev et al.l 
dl982h . the emission height of various radio frequencies in the 
magnetosphere of the pulsar can be derived using the frequency 
dependence of the component separation. Our analysis shows 
that a y-°- 30±002 relation is valid for the component separation in 
the frequency range considered here and is s lightly steeper than 
the y~° 24 and y~°- 2 6 scaling laws derived by iBackerl d 19721) and 
Siebe r et all d 19751) . respectively. 
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i1 1 12 + 50 



11 1 33+1 6 



B0031 -07 



2.05 




Pulse Phase 

Fig. 1. Average pulse profiles of the three pulsars. The profile 
changes with sky frequency is clearly evident. The profile for 
PSR Bl 133+16 at 163.5 MHz is not shown because it was too 
corrupted by RFI. The profiles are shown at 809, 580, 460 fis 
time resolution for PSRs Bl 112+50, Bl 133+16 and B0031-07 
respectively. All profiles are normalised by the maximum flux 
density in the band and are aligned at the phase of maximum 
intensity in each band. Profiles at 116.75 and 130 MHz have 
comparatively more noise because up to ~15% of pulses were 
removed to reject RFI. 



3.4. Radio spectrum 

The observations presented here span a reasonably wide fre- 
quency range and we use the flux of the pulsar in each of the 
recorded bands to compute the pulsar spectra. The pulsar flux 
used here is the mean continuum flux density of these sources, 
computed from the flux-calibrated pulse profiles. The functional 
form of the radio spectrum of the pulsars was assumed to follow 
S(v) oc y", where S(v) is the pulsar flux at frequency v and a is 
the spectral index. The parameter a was estimated using a least 
squares fitting procedure and the result is displayed in Fig.[2]and 
Table [2] for the three pulsars. To account for the effect of RFI 
and possible errors in telescope parameters, a 10% uncertainty 
in the estimation of the pulsar flux density in each band is con- 
sidered, which dominates other errors in S/N estimates because 
all pulsars are strongly detected. 



in 
6 



CN 

d 



CD 

Q 



LO 

d 



in 
d 




2.05 2.1 2.15 2.2 

Log Frequency (MHz) 



2.25 



Fig. 2. Spectra for three pulsars in the frequency range 100-200 
MHz. Errors bars correspond to a 10% uncertainty in the pulsar 
flux estimation. The best fit power law curves are also displayed. 
For PSR Bl 133+16, a broken power-law provides a closer fit to 
the measured fluxes. The flux in the RFI-affected band centred 
at 163.5 MHz is not included for PSR Bl 133+16. 



We now consider the effects of scintillation on the estima- 
tion of the spectral indices. The length of these observations are 
between 5 and 9 hours (see Table. [TJ and may be susceptible to 
scintillation effects on these timescales. Diffractive scintillation 
does not influence our results here as we now demonstrate. The 
frequency-phase plots displayed in Fig. [3] shows that the inten- 
sity modulation induced by this type of scintillation is confined 
to a few channels. To help validate this, the diffractive scintilla- 
tion bandwidth and time were derived from values published in 
other work (see Table. [2), using Avdiss 00 y4 ' 4 an d ^toiss K v 12 
scaling laws. From the foregoing discussion, it is clear that the 
observed and derived diffractive scintillation bandwidth is nar- 
row for these pulsars at low frequencies and that we average 
over many "scintiles" with the wide band used here. The ef- 
fects of the refractive scintillation should also be considered for 
th ese long observations, and this timescale is derived as outlined 
in [Stinebring & Condon! d 19901) and listed in Table. [2] Hence, 
it is clear that the diffractive scintillation does not affect the 
results derived here. The refractive scintillation may affect the 
these long observations, but are accounted for in the 10% uncer- 
tainty in flux estimation considered here. The derived timescale 
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Table 2. Measured and derived parameters of three pulsars. 



Pulsar 


S " 




»2 


Avdiss c 


TV- 


77 


Reference 




mjy 






kHz 


sees. 


days 




B1112+50 


460 + 46 


1.12+0.56 




77 - 447 


115-187 


2.6-1.1 


1 


B1133+16* 


1030 + 51 


2.33±2.55 


-3.81+2.24 


7-40 


48-77 


11.8-4.9 


2 


B0031-07 


821 + 41 


0.71+0.45 




3-19 


212 - 343 


109.1 -45.2 


3 



Notes. 

<a> The flux density quoted here is the value measured at 142.25 MHz. 

( *' A broken power-law spectrum results in a bett er fit for this pulsar. 

(c) Scintillation parameters are derived from - (D lCordes & Rickettll 19981 : (2l lKramer e"t"alll2003r . Gl lStinebring & Condonll 19901 



T r shows that refractive scintillation, if there is any, affects only 
the spectral index estimation for PSR Bl 1 12+50. 

PSR Bl 1 12+50 shows a slightly positive spectral index and 
we find a value of a = 1.12 ± 0.46. Even though this pulsar is 
the weaker of the three pulsars, it is detected with high signifi- 
cance, as is clearly evident from Fig. Q] The somewhat poor fit 
to the pulsar spectra in these observations is not caused by fluxes 
corrugated by diffractive scintillation for the reasons explained 
above and may have its origins in the residual RFI and/or re- 
fractive scintillation. The results presented here imply that the 
spectra turns over at a frequency higher than 180 MHz. 

For PSR Bl 13 3+16, the middle panel of Fig.[2]can be com- 
pared to Fig. 4e in Deshpande & Radhakrishnan] d 1992b . which 
extends down to 34.5 MHz. Their figure can be interpreted as 
indicating a spectral turnover in our frequency range. With this 
consideration, the spectrum was modelled by a two component 
power-law defined as 



S(v) 



C 2 v° 



if v < Vb 
if v > v/,, 



(2) 



where, Vb is the break frequency, C\ is a constant, and normal- 
isation yields C 2 = Civ^ 1- " 2 . The values of break frequency 
of 139.75 MHz with spectral indices a\ = 2.33 + 2.55 and 
a2 = -3.81 + 2.24 were derived from the best fits to the data. 
The effect of scintillation gives rise to "scintiles" with typical 
widths of < 40 KHz, which tend to average out within each 2.5 
MHz band and do not contribute to the intensity variations in the 
1 10-180 MHz frequenc y range. Co mbined with the 75% error in 
400 mjy at 151 MHz dSieberlll973h and 900 mjy at 1 1 1 MHz 
(Figure 3 in iMalofeev et alj 19941) . we conclude that our results 
are consistent with those reported previously. 

We find a relatively flat spectrum for PSR B003 1 -07 with an 
index of a = 0.71+0.45. As for the other two pulsars, B003 1 -07 
has partially resolved scintillation (see Fig. [3]) at these frequen- 
cies that do not con tribute to the es t imated fluxes. For sky fre- 
quencies v < 1 GHz. IMalofeev et all d 19941) find that the spectral 
index of this pulsar is a = —1.2 ± 0.26. Consideri ng the two data 
points below v = 100 MHz and the error limits in Mal ofeev et alJ 
(1994), the spectral index derived here is consistent with their 
work. The trend seen in the spectra of this pulsar implies that the 
spectra turns over at a frequency higher than 180 MHz. 



4. Bright single pulses 

Several bright pulses were detected from all three pulsars and 
examples are shown in Fig. [3] These pulses are very narrow and 
often less than 1 ms in width. In general, giant pulses are char- 
acterised by their widths, the phase of the pulses with respect 
to the average emission profile, and the pulse energy (Knight 



. The classical giant pulses are al so found to be broad band, 
for example in B 1937+21 Av/v ~ 0.2 dPopov & Stappersl2003l) . 
where v and Av are the observation frequency and bandwidth. In 
this study, we find similar bright, broadband bursts from the pul- 
sars studied, although they are not as energetic as the classical 
giant pulse emitters. 



Table 3. Brightness temperature of single pulses. 



Pulsar 


d 


T B 




kpc 


K 


Bl 112+50 


0.32° 


3.51 x 10 26 


B1133+16 


1.06 6 


1.93 x 10 27 


B0031-07 


0.35 r 


1.32 x 10 27 



Notes. 

<n) DM distance from lTavlor & Cordesll 19931. 
® Parallax distance from|Chatterjee_et^22009. 
w Parallax distance from Brisken et al. 2002. 



The implied brightness tem perature of the single pulses 
can b e computed from (pp. 79; lLorimer. D. R. and Kramer. M.I 
2005), 



T R = 



S peak \l vAt 

2nkB j\ d 



-2 



(3) 



where S pea k is the peak pulse flux, ks is the Boltzmann's con- 
stant, d is the distance to the pulsar, v the sky frequency and At 
is the pulse width. Using At values from Tableland the known 
distances for these pulsars, the computed brightness tempera- 
tures are listed in Table [3] 

While the brightness temperatures of these pulses are sev- 
eral orders of magnitude lower than the 10 35 - 10 40 K seen in 
typical giant pulses from the Crab pulsar or PSR B9 137+21, 
the large values of Tg do point to a coherent emission mecha- 
nism for these pulses. From the frequency-phase plots in Fig. [3] 
these bright pulses clearly represent broadband emission, with 
Av/v = 0.3, and can be compared to the value of 0. 1 and 0.2 for 
the Crab pulsar and PSR B 1937 +21 at 1400 MHz. The scintiles 
in all three pulsars are also visible in this figure, though they 
are more evident in PSRs B1133+16 and B0031-07. Therefore 
these bright pulses are not due to the occasional focussing ef- 
fects brought about by scintillation. The similarities between the 
bright pulses found here and the classical giant pulses are indica- 
tive of high energy counterparts as in the Crab pulsar and PSR 
B1937+21. 
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B1 1 1 2 + 50 B1133+16 B0031-07 




Pulse phase 

Fig. 3. Narrow bright pulses from the three pulsars detected at all eight sky frequencies. The total intensity displayed in the top 
panels is computed by summing the signal in the 8 bands. The lower panels show a dedispersed pulse as a frequency-phase plot. 
The abscissa in these plots correspond to the pulse phase at 116.75 MHz. Some residual RFI is visible as low-level intensity 
modulating the entire pulse phase in the frequency-phase plots of all three pulsars. 



5. Dispersion measure variations 

Pulsar dispersion measures are normally estimated 
from multi-frequenc y mu lt i-epoch timing observations 
dPhillips & Wolszczanl Il992t iBacker et all 1 1993b . However, 
in some pulsars the absence of an appropriate fiducial point in 
the average emission pro f ile int roduces errors in the estimated 
DM values. lAhuia et alJ 12001 used cross-correlation of the 
pulse profiles at two different frequencies of observation to 
partially overcome this problem. These authors also point 
out that an alternative means of estimating the DM is to use 
multi-frequency simultaneous observation of the pulsar, and 
use a fiducial point in the average emission profile to align the 
profiles at the two frequencies. Our detection of single, narrow, 
and bright pulses at low sky frequencies provides in an excellent 
reference to align the pulses and hence compute the DM very 
accurately for these low DM pulsars. The frequency resolution 
of 64 channels across each 2.5 MHz band gives a total of 512 
channels in the 1 10-180 MHz range, which permits the accurate 
determination of the quadratic dispersion curve. We generally 
started with the catalog DM, which was varied until a peak S/N 
in the combined pulse profile was attained. The frequency-phase 
plots were then visually inspected for the pulse alignment using 
the new value of DM. To compute the change in the DM at the 
epoch of our observations, the dispersion smearing caused by 
the use of an incorrect DM within a single frequency channel 
was considered negligible. 



From their multi-year timing data, iHobbs et al.l (120041) 
find DM=4.864 cm _3 pc for PSR B1133+16 corresponding 
to the epoch MJD 46407 and also determine d(DM)/df « 
-0.0008 cm^pcyr -1 . Our analysis shows that the DM has 
changed to 4.844+0.002 cirT 3 pc at epoch MJD 54793 from our 
reference value of 4.871 cirT 3 pc at epoch MJD 41665. The er- 
ror estimate comes from the DM uncertainty caused by the mis- 
alignment of a broadband pulse by one phase bin across the 
110-180 MHz band. The total DM c hange, ADM , found here 
is consistent with the predictions of IHobbs et al.l (|2004). The 
somewhat high yearly rate of change in the DM of this pulsar 
probably arises from its large proper motion ~630 km-s -1 . The 
ADM= 0.0265 cirT 3 pc corresponds to a 0.5% change in DM 
and failing to correct for this change leads to the peak flux of the 
bright pulses being underestimated by ~3%. 

For PSR Bl 112+50, we find that the DM increased from 
9.16 cnr 3 pc at an epoch MJD 50899 to 9. 1750+0.0001 cm 3 pc 
at epoch of MJD 54795 resulting in ADM=0.015 cm 3 pc. The 
error in our DM value is derived from the pulse misalignment by 
one p hase bin in a 2.5 M Hz band at 147.5 MHz. We again re- 
fer to lHobbs et al.1 (12004) for DM estimates of this pulsar, where 
they find DM=9.195cnT 3 pc and a yearly rate of d(DM)/df = 
-0.004 cirr 3 pcyr~'. These values are inconsistent with the rate 
of -0.0014cirT 3 pc yr~' derived from our analysis. One reason 
for this discrepancy could be the difference between the high 
frequency observat ions at 408-1600 MHz used in the work of 
Hobb s et all (|2004) and the much lower frequencies used in our 
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study which are more sensitive to DM changes. The second rea- 
son could be that a single value of d(DM)/df is insufficient im- 
plying that higher order terms of the DM time derivative may be 
necessary to accurately model the DM variations. 

We do not detect a change in the DM of PSR B0031-07 
from the reference DM value of 10.89 crrT 3 pc. The narrow and 
bright single pulses have no discernible errors when the sig- 
nal is dedispersed with this DM. However, the pulsar's DM of 
10.89 cm~ 3 pc is inconsistent with the value of 11.29 cirT 3 pc 
at this epoch expecte d from the -0.007 cirT 3 p c yr' 1 listed 
et alj d2004l). Mo r eover, in th e studies of iTaylor et alj 
uzmin & Ershovl ([2004), and ISmits et alJ (l2007h . the 



authors use a DM of 10.89 cm pc, which implies that the DM 
of the pulsar has not changed since 1993. Therefore both the DM 
value and the rate of DM change for this pulsar in Hobbs ~et~aT] 
(2004) may be incorrect. 



6. Single pulse analyses 



Our long observations resulted in the detection of a large 
number of pulses and are shown in Table Q] Robust statistics 
can be obtained by computing the distribution of single pulse 
energies E and peak intensities S p . A scatter plot of the peak 
fluxes and the average emission profile for the three pulsars are 
shown in Fig. |4] The 2.5 MHz band at 156 MHz was chosen out 
of the eight bands for further study because it showed compara- 
tively less RFI. We note that after an exhaustive RFI search and 
removal, all eight bands show similar trends. 

The presence of narrow bright pulses was evident when the 
data was examined visually in the RFI cleaning stage. These 
unusual emission entities should be evident in the single pulse 
energy and flux distributions. For example, giant pulses in the 
emission of the Crab pu lsar should result in a power law in the 
pulse flux distribution (Lundgren et al.lll995l) . For the three pul- 
sars considered here, the distribution of the pulse flux and energy 
are shown in Figs. [5] and [6] respectively. The pulse energy was 
computed by summing the emission in a window of width equal 
to the average emission profile. The single-pulse peak flux was 
computed using the product of the peak S /N and S av (Eq.[TJ. 
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Fig. 4. Single bright pulses and the integrated pulse profiles at 
156 MHz for the three pulsars. The flux scale on the left ordinate 
corresponds to the peak flux of the single pulses and the scale on 
the right ordinate axis corresponds to the average pulse profile. 
The phase displayed here is different from that of Fig.[3]due to 
the signal being dedispered to a different central frequency. Peak 
flux densities of the single pulses are shown as crosses. Only 
pulses > 20cr are shown in these plots which represent a low 
flux cut-off in the fluxes of single pulses. 
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Average Pulse Energy ( E/<E> ) 

Fig. 5. Pulse energy distribution for the three pulsars. The ab- 
scissa shows pulse energy normalised to the mean pulse energy. 
Dark and light grey lines correspond to the on-pulse and off- 
pulse energy distributions. The clear peak at zero energy is in- 
troduced by the RFI affected pulses replaced by zeros. 
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Fig. 6. Plot of the peak flux distribution of the single pulses for 
the three pulsars. The dashed line shows 10 x S aViPea k for each 
pulsar. 



6.1. PSRB1112+50 

The top panel of Fig.|4]shows that this pulsar emits bright pulses. 
Some of these pulses show considerable structure, with narrow 
subpulses separated by a few milliseconds. These pulses have a 
slight preference to the trailing phase of the average pulse pro- 
file, although a clear phase relationship is not evident. This is 
simil ar to the giant pu l se-lik e emission from this pulsar reported 
bv lErshov & Kuzminl(l2003l) . Fig.[5]shows the pulse energy dis- 
tribution and the bright pulses with energies greater than IQx(E) 
comprise the trailing part of the distribution. An indication of the 
null emission in this pulsar is visible from the slight turn-over in 
the pulse energy distribution near 3 x (E), which is similar to 
the offpulse energy distribution. Our data is consistent w ith the 
> 60% nul ling fraction at 1 412 MHz dWrightet al.ll 19861) and at 
408 MHz (Ritc hingslll976h . If the null emissions are accounted 
for with this nulling fraction, then (E) increases by 60%. This 
causes a decrease in the energy of bright pulses in the 10-13 
x(E) range. The pulse flux distribution also displays a shallow 
power-law like tail in the peak flux distribution (top panel in Fig. 
|6| produced by the occasional bright pulses in this pulsar. 

6.2. PSRB1133+16 

iKramer et al.l d2003l) reported bright pulses from this pulsar at 
a frequency of 5 GHz but not at 1.4 GHz. We detected bright, 



narrow pulses in all eight bands in the range 1 10-180 MHz with 
peak fluxes greater than 10 times the peak of average flux density 
as seen in the middle panel of Fig. [6] The pulse energy distribu- 
tion shown in the middle panel of Fig. [5] can be qualitatively 
approximated by a power law. The distribution does not exhibit 
pulses with energies >10 x (£), which is the normal working 
definition of giant pulses (Cogna rd et al.l[l996h . From the mid- 
dle panel in Fig. [4] a preference for the bright pulses to fall 
on the inner edge of the two components is evident. Some of 
these bright pulses contain two distinct and narrow subpulses. 
Brigh t pulses from this p ulsar at 327 MHz were also noted by 
Herfindal & Ranki n| d2007). The turn-over in the pulse energy 
distribution ~1.0x(.E) is indicative of abimodal distribution, and 
hence of nulling. The pulsar show s a ~20% null fraction at 327 
MHz (Herfindal & Ran idnll2007l) . but at our frequency we are 
unable to identify the null distribution because we were limited 
by the system sensitivity. 

6.3. PSRB0031-07 

iKuzmin & Ershovl (|2004|) reported giant pulses from PSR 
B0031-07. Our observations reveal that the bright and spiky 
pulses do not satisfy the 10x(£') criterion but are easily greater 
than 10 x S aVyPea k, as seen in Figs. [4] and [6] Fig. |4] also reveals 
the emission of bright puls es on the leading edge of t he average 
pulse profile. The work of lErshov & Kuzm in (2005) (see their 
Table 1) suggests that the bright pulses in this pulsar have a steep 
spectral index. We also d raw the reader's attention to Fig. 4 of 
IKuzmin & Ershovl d2004l) . where the authours report the detec- 
tion of "double giants" in this pulsar at 40 and 111 MHz. These 
pulses were detected in our observations, as seen in the rightmost 
panel of Fig. [3] These bright subpulses with two distinct peaks 
show structure spanning a few milliseconds implying that they 
are not microstmcture in the emission. As with PSR Bl 133+16, 
a turn-over in the pulsar energy distribution ~1.5x(£') is indica- 
tive of null emission. 

7. Microstructure 

Several single pulses from both PSR Bl 133+16 and B0031-07 
show pulse structures considerably narrower than the width of 
the average pulsed emission. The presence of pulses with mi- 
crostructure became evident when several thousand pulses were 
examined in the process of cleaning RFI from the data. To 
quantify the micros tructure, the autocorrelation function (ACF) 
dCordes et alj [l990) of the single pulses and their averages were 
computed. The discrete ACF is defined by 

nbin 

R(T) = J]l(f)-I(t + T), (4) 

t=0 

where I{t) is the intensity of the pulsar. Fig. [7] shows R(j) for a 
single pulse and the averaged ACF of several pulses. 

For PSR Bl 1 12+50, the single pulses were quite narrow, de- 
void of any feature that can be recognised as microstructure. 
Therefore, this pulsar is not considered any further in this sec- 
tion. The separation between the two main emission components 
can be clearly seen for Bl 133+16 and is the feature at correlation 
lag ~30 ms. The ACF for the pulse shown in the uppermost panel 
exhibits the periodicity of microstructure (i.e., the local maxima 
near correlation lag 5 ms in the middle panel of Fig. [7]). It is 
known that PSR Bl 133+16 shows structure down to ~8/is and a 
characteristic intensity fluctuation of the order of ~575 - 663 /is 
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50 
time lag (ms) 

Fig. 7. Intensity autocorrelation functions (ACF) of pulses with 
structure for two pulsars. The uppermost panels present exam- 
ple pulses. The autocorrelations displayed in the middle panel 
correspond to the pulses above them. The two bottom panels 
display the averaged ACFs for 120 and 132 pulses for PSRs 
B0031-07 and Bl 133+16 respectively. The averaged ACF's for 
two bands centred on 116.75 MHz and 173.25 MHz (dashed 
line) are shown. 



in the 1 1 1-196 MHz range dHankinsll 19721) . Our final resolution 
of 580 fj.s was sufficient neither to resolve narrower structures 
nor to permit us to estimate the persistent quasi-periodic inten- 
sity variations in the averaged ACF. 

In the case of PSR B0031-07, single pulses display nar- 
row bursty features. However, the single pulse ACF in the right 
middle panel of Fig. [7] shows no characteristic variations im- 
plying that the pulse consists of non-periodic microstructure. 
These nonperiodic features were also found in PSR B0950+08 
dLange et alj fl998). The average subpulse separation is seen in 
the averaged ACFs and is visible at -55 ms for PSR B0031-07. 
Expressed in pulse longitude, this is 20°, which is also the value 
of P2 found in Sect. [8] For the average ACF's displayed in the 
two lowest panels of Fig. [7] the pulses chosen are detected si- 
multaneously in the two bands. The slight shift of local maxima 
to the left cannot be measured with high significance, implying 
that the subpulse separation decreases with frequency as can be 
se en by comparing t h e subpu lse p roperties at 328 and 1380 MHz 
in IWeltevrede elatl d2006al) and IWeltevrede etail (120071) . The 
slight break at -10 ms in the ACF of the single pulse represents 
the width of microstructure in the pulsar. The reason for the lack 
of a prominent characteristic timescale in the microstructure of 
PSR B003 1 -07 may be our coarse time resolution. 



8. Fluctuation spectra 

Since the three pulsars under discussion display bright and 
bursty emission, it is natural to expect considerable pulse- 
to-pulse variation. In some pulsars (eg. PSRs B0809+74, 
B 1 944 + 1 7), drifting and per sistent microstructure were detected 
leading Corde s~et alj d!990t) to suggest that these two phenom- 
ena are correlated. Moreover, PSRs Bl 133+16 and B0031-07 
show drifting subpulses at other frequencies. To ascertain the 
drifting behaviour at these frequencies, we undertook sub-pulse 
drift analysis. 

To perform the analysis presented in this section, the 
pulse stacks were subjected to a pulse drift analysis 
based on the longitude-resolved modulation index (LRMI), 
longitude-resolved fluctuation spectrum (LR FS), and two- 
dimensional fluctuation spectru m (2DFS) dEdwards & Stapperi 
I2002t IWeltevrede et al.l l2006al) . The sliding two-dimensional 
fluctuation spectrum (S2DFS) was also computed to identify the 
onset of the various drift modes (Servla k et al.l2009T) . The result- 
ing plots from the LRMI, LRFS, and 2DFS analysis are shown 
in Fig. [8] In this plot, only LRMI values that are detected with 
a significance of more than 6cr are shown (top panel; Fig. [SJ. 
In the LRFS and 2DFS plots, drifting subpulses produce regions 
of enhanced brightness called "features". The vertical drift rate 
P3 and the horizontal separation P2 of the subpulses are defined 
to be the centroid of a rectangular region in the 2DFS contain- 
ing the features (Eq. 6: IWeltevrede et alj |2006a). The sign of P2 
denotes the direction of the drift, i.e., whether subpulses drift to- 
wards the leading (negative) or trailing (positive) edge of the pul- 
sar's average profile, respectively. In what follows, the S2DFS 
analysis results are discussed only for PSR B003 1 -07. The anal- 
ysis of PSRs Bl 133+16 and Bl 112+50 using the S2DFS tech- 
nique provided results that are comparable to those obtained 
from the 2DFS method alone. 

The results presented in this section are for the 2.5 MHz 
band centred on 156 MHz, as in the previous section. In the RFI 
cleaning procedure discussed in Sect. 13.11 fewer than 2% of the 
pulses were replaced by zero, and therefore the effect is consid- 
ered small enough to be neglected. The drift analysis was per- 
formed on all eight bands, and the results presented here are rep- 
resentative of the other bands. Within the 1 10-180 MHz range, 
there are no noteworthy changes in drift features. The modula- 
tion index was also in general agreement with the results found 
in the band reported here but are higher than the values reported 
in previous high frequency studies (see Table[4]l. 



8.1. PSRB1112+50 

This pulsar shows a broad, non-specifi c vertical drift featur e as 
seen in the panels on the left in Fig. [8] Wright et al. (1986|) re- 
ported the presence of three different drift modes in this pulsar 
and that the dominant mode in the trailing part of the average 
pulse profile showed a periodicity of P3 ~ 6Po- At 328 MHz, 
Welt evrede et al.l (I2007I) report that Po = 40~Io degrees and a 
P 3 = 9 + 5P . We detect P 2 = 5^ de grees and subpulse modula- 
tion of P3 = 6 + 3fo- We also note a greater amount of power in 
the entire range of P3 centred around P0/P2 = 0. While the P3 
found here is consistent with values reported elsewhere within 
the errors, the di s crepan cy between the P2 results reported by 
Weltevrede et al. I J20071) and our own can be explained by the 
higher S/N of individual pulses in our observations. This results 
in a more prominent feature in the spectrum and less convolution 
with the higher power in the entire range of P3 centered around 
P0/P2 = 0. 
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Fluctuotion frequency (P /P 2 ) 

Fig. 8. Plots of the fluctuation analysis for the PSRs B1112+50, B1133 + 16, and B0031-07. The average pulse profile, and the 
longitude-resolved modulation index (LRMI) are shown in the uppermost panels. The panels directly below the average pulse 
profile present the longitude-resolved fluctuation spectrum (LRFS). The two-dimensional fluctuation spectra (2DFS) are shown in 
the third and fourth rows. For PSR Bl 133+16, the 2DFS plots are shown separately corresponding to the two components of the 
average pulse profile. The panels on the left and bottom of the LRFS and 2DFS images contain the vertically and horizontally 
integrated spectra, respectively. 



It is possible that only the prominent drift mode active at 
1412 MHz (mode 2 of Wrig ht et allll986l) is also active at this 
frequency. The support for this comes from the similarity of 
the average profile at our frequency to the average profile com- 
puted from 193 pulses containing drift mode 2 in [Wright et al.l 
d!986l) . In their work, the profile computed from 176 pulses 
in mode 1 displays a prominent double peak. From the LRMI 
plot, the modulation index m ~ 3.5 is quite high at this fre- 
quency. Combined with the values reported elsewhere (see Table 
@J, this indicates that the pulse-to-pulse intensity variations in- 
crease with decreasing frequency. This again suggests that the 



large amplitude pulses are more common in this pulsar at lower 
sky frequencies. The LRFS shows an excess near P0/P3 ~ 
owing to th e large nulling f raction of ~60% in the emission of 
this pulsar (lRitchingslll976l) . which can be also seen in Fig. [5] 
However, this can only be indicative of nulling because some 
contribution to the excess seen comes from the zeros introduced 
by the RFI cleaning procedure. Therefore, these plots cannot be 
used to measure the pulsar's nulling fraction. 
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Table 4. Drift analysis results for three pulsars. 



Pulsar 


P 3 


Pi 




77! 






[Po] 


[deg] 


21cm 


92cm 


192cm 


B1112+50 ... 


6±3 


5 +4 

-2 


1.5 


2.1 


3.5 


Bl 133+16° . . . 


30 ±5 
34+5 




1.4 


0.8 


2.2 


B0031-07* . . . 


12 2 +0 - 2 
6.7 ±0.1 


-21!? 

-21 6 +a5 

z ' 1- -0.2 


1.2 


1.4 


3.1 



Notes. Vertical and horizontal separations of the subpulses are listed in 
Cols. 2 and 3, respectively. The value of modulation index m is the mini- 
mum in the longitude-resolved modulation index and is shown for three 
differ ent wavelengths in the las t three columns. The v alues in Col. 4 are 
from [y Velte vrede et al. (2006a), values in Col. 5 from Weltevred e et alJ 
j2007fK and the results in the last column are from this work. 
(o) Values presented for the leading and trailing component. 
(b) Values presented for the drift mode "A" and "B". 



8.2. PSRB1133+16 

The results of the analysis presented here confirms the presence 
of the long-period feature P3 ~ 30 at this frequency, which 
was previously reported at frequencies greate r than 300 MHz 
(IWeltevrede etal] [20061 iHerfindal & Raiikinll2007h . In Fig. [8] 
(middle panel), the horizontal component of the drift feature, P2 
is very weak and consistent with no ordered horizontal separa- 
tion of the subpulses. The feature near P0/P3 = 0.18 is present 
in both components of the pulse profile and the LRMI plot shows 
the presence of a significant modulation in the bridge emission 
region. We find P3 = 30 + 5Po for the leading component and 
P3 = 34 + ^Po for the trailing component. Apart from the afore- 
mentioned P3 feature, the LRFS shows a non-ordered feature 
that modulates the emission throughout the entire stretch of the 
analysed pulse sequence. This is also evident from the elevated 
baseline of the collapsed spectra shown in the side panels of Fig. 
[8] The excess power seen near P0/P3 = is due to the combi- 
nation of null emission and the pulses replaced by zeros in the 
process of cleaning RFI, though as in PSR Bl 1 12+50 this cannot 
be measured reliably from this analysis. 

8.3. PSR B0031-07 

Th e three different drift m odes in this pulsar were first reported 
by iHuguenin et alJ (ll970l) of which two are clearly visible in 
our analysis (Fig. [8] right-hand panel). The very low frequency 
feature at P0/P3 ~0.02 comes from the nulls present in the 
pulse stack. The nulls are a combination of the null emission 
and the zeros introduced by RFI cleaning procedure. We detect 
this feature in all the eight observed bands. The feature seen 
at P0/P3 = 0.15 denotes the prominent drift mode "B" that 
introduces harmonics and hence appears as the feature seen at 
P /p 3 = 0.3. The other feature seen at P /P 3 = 0.08 is the drift 
mode "A", which can be identified with regular upward drift- 
ing bands in the pulse stack. The drift mode "C" is not seen 
in the 2DFS spectrum. We used the S2DFS method to analyze 
the occurrence rate of different drift modes. While modes "A" 
and "B" were clearly detected in the S2DFS output, mode "C" 
remained undetected. The mode "C" has a low occurrence rate 
and is active only for ~2.6% of the obs ervation tim e at 32 7 MHz 
dVivekanand & Joshill 19971) . However. Isinits et all (120051) detect 



this mode marginally at 328 MHz but not at 4.85 GHz, indicat- 
ing a steeper spectral index for the drift mode. Therefore, other 
than susceptibility to RFI and lower sensitivity, it is unclear why 
this mode is not detected at these frequencies. 

9. Discussion 

The long observations reported in this paper provide good sen- 
sitivity to persistent weak features in the average emission of 
the pulsars. We were thus able to test the hypoth esis of in- 
terpuls e emission in PSR Bl 1 12+50 proposed by Wright et al. 
(1986) at 1400 MHz. Based on the outer gap emission model 
then proposed, these authors recommend the use of observa- 
tions at low sky frequencies to detect the inter pulse, as the 
model predicts large fanbeams. Our deep exposures at low sky 
frequencies shows no evidence of inter-pulse emission with an 
upper limit of 3.5 mJy at 156 MHz (see Fig. [TJ, which is < 
0.6% of the main pulse peak flux. The non-detection of an in- 
terpulse in this pulsar is unsurprising given that only ~1.5% of 
the 184 7 known pulsars possess an i nterpulse emission com- 
ponent (IWeltevrede & Johnstonll20 08i). Examination of polari- 
metric profiles at 400 and 608 MHz from the Jodrell Bank 
Observatory archives shows that the position angle swing is too 
shallow, which would result in a large impact angle to the ob- 
served radio beam. Our observations and the indication from 
past polarisation studies shows that if the neutron star has a dipo- 
lar magnetic field, PSR Bl 1 12+50 is not an orthogonal rotator. 
Consequently, the alternate magnetic pole never sweeps across 
our line of sight. 

The pulsar spectral indi ces found in this study indicate a 
turn-over near 180 MHz. iMichell d 19781) modelled the pulsar 
emission as consisting of a coherent part and an incoherent syn- 
chrotron component. This model also predicts a spectral index of 
-8/3 and a spectral turn-over below 100 MHz based on the finite 
energy of the injected particles, and that the radiation is domi- 
nated by incoherent emission mechanisms. However, this con- 
trasts with the presence of bright pulses and clear single pulses 
near the putative turnover frequency implying th at the radia- 
tion is still of a coherent nature. The simulations of lGil & K iiak 
(119921) . which assume that the pulsar radiation is produced by 
emission cones, clearly illustrate that the estimated intensity is 
very sensitive to the viewing geometry. If the radiation originates 
quite high in the magnetosphere (as at low sky frequencies), it is 
possible that the viewing geometry changes quite rapidly. This 
in turn can give rise to the observed pulse flux near the spec- 
tral turn-over resulting in spectra that does not vary smoothly. 
Detailed polarisation studies at these frequencies can be used 
to test whether the observed frequency dependence of the pulse 
flux has its origins in the viewing geometry, although the pulse 
shape is quite stable. In summary, past studies have been in- 
conclusive on the value of spectral index for PSRs Bl 133+16 
and B0031-07 at low sky frequencies. The results presented 
in this work suggest a spectral turn-over near ~140 MHz for 
PSR Bl 133+16 and above 180 MHz for PSRs B0031-07 and 
Bl 112+50. 

As pointed out by other authors, the discrepancies in pulsar 
DMs stem from the difficulty in aligning multi-frequency pulse 
profiles and the choice of a fiducial p oint in the average profile 
(iHobbs et alJf2004HAhuia et~aj]l2007l) . The narrow bright pulses 
at low sky frequencies in the pulsars studied in this work have 
permitted a very accurate estimation of the DM. We note that 
scatter broadening of the emitted pulses can limit the accuracy 
of the DMs estimated by both other methods and ourselves. This 
can be another source of the DM discrepancies. Because of the 
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broad band nature of the pulsar emission mechanism, we cannot 
rule out the claim that some pulsars may have slightly differ- 
ent values of DM in d ifferent parts of the frequency spectrum 
(Kardashe v et al.lfl982h . However, with the S/N in this work we 
do not measure any deviation from the v 2 dependence of the 
interstellar medium contribution to dispersion in the 110-180 
MHz band. Simultaneous broad-band observations with the ex- 
isting instruments and LOFAR can address these outstanding 
questions effectively. 

The analyses of the single pulse energy distributions revealed 
that only PSR Bl 1 12+50 exhibits pulses of energy > 10 x (E). 
The presence of null emissions in these pulsars leads to an under- 
estimate of the average pulse energy. Therefore, this is an impor- 
tant consideration when interpreting pulse energy distributions 
for pulsars that exhibit the nulling phenomena. The width and in- 
tensity of the bright pulses in this study are quite similar to those 
of the normal giant pulses suggesting that they are caused by 
similar emission mechanisms. Sources of classical giant pulses 
display a phase correlation between the giant pulses and the 
high energy non-ther mal X-ray emissio n, as in the Crab pulsa r 
dLundgrenetal.il 19951) . PSR B 1937+21 ( Cusum anoetaT1l2003h . 
and in other millisecond pulsars (iKnight et a l. 2006a). However, 
these are energetic pulsars with large values of Bic- It is known 
that P SR Bl 133+16 shows X-ray emission (Kargaltsev et alJ 
2006), which the authors claim contains both thermal and non- 
thermal components. With their coarse time resolution of 3.2 
s and low photon counts, these authors were unable to detect 
any X-ray pulsations. However, the non-thermal high energy 
emission may be related to the bright pulses observed in our 
work. For PSRs Bl 112+50 and B0031-07, a deep search in 
the archival X-Ray data might reveal any correlation that might 
exist. The bright d ouble -pulses reported in this work and by 
iKuzmin & Ershovl (120041) are indicative of a global magneto- 
spheric phenomena, despite the different physical locations of 
the emitting regions. However, it could also be that drifting oc- 
casionally brings these two regions into our line of sight result- 
ing in the double pulses. Detailed analysis of these double pulses 
will provide insight into the phenomena that might be active over 
the whole magnetosphere. 

The preference of the bright pulse emission to occur at the 
phase coincident with the peak in the average emission profile 
is a feature common to all pulsars in this study. This argues that 
these pulses are preferentially emitted further away from the last 
open magnetic line in the ma gnetosphere, in contrast to pul - 
sars such as PSR B 1937+21 dKinkhabwala & Thorsed 12000). 
The bright pulses in B0031-07 are more frequent and bright 
at even lower sky frequencies when we com pare our results in 
the 11 0-180 MHz range with those found by Ershov & Kuzmin 
(2005) at 40 MHz. This is also valid for the other two pulsars 
studied here - bright pulses are more common in these pulsars 
at lower frequ encies. This finding contrasts with the models of 
iPetroval d2004l) in which the bright and narrow pulses at high fre- 
quencies were explained by inverse Compton scattering of low 
frequency radio photons. Simultaneous multi-wavelength radio 
observations might help to resolve this issue. 

In the pulsar models proposed by iRuderman & Sutherland! 
(fl97l . the radio waves are emitted tangentially to the diverging 
magnetic field lines by the relativistic electron-positron pairs. 
In this model, the width of the observed average profile is then 
defined by the last open magnetic field lines. Furthermore, the 
model predicts that the emitted radio frequencies v follows a 
v oc r~ 3/2 relation, where r is the radial distance from the cen- 
tre of the star. This leads to the lower frequencies being emit- 
ted higher in the magnetosphere than the high frequency ra- 



dio waves, giving rise to a radius-to-frequency mapping (RFM). 
Our observations show a decrease in the component separa- 
tion in the average pulse profile of PSR Bl 133+16 with fre- 
quency, which scales as v~ 3 . After correcting for the residual 
DM smearing, the retardation and abberation time delays pre- 
dicted by the RFM model can be measured from the 116 and 
173 M Hz average profiles. Using a method similar to Cordes 
(fl978l we find a lower limit of 5.6 x 10 5 m for the height 
of 116 MHz radio emission. Depending on the methods used, 
previously measured emission heights were < 10 5 -10 5 8 m and 
< IP 4 7 - 10 5 - 8 m for the radio emission in the 40- 1 400 MHz range 
(ICordesl[T97l iMatese & Whitmird [19801) . The upper limit in 
these works corresponds to the 40 MHz radiation emitted at ra- 
dius r < 10 5 8 = 6.3 x 10 5 m. Therefore, the derived emission 
height of 5.6 x 10 5 m for the 116 MHz radiation is consistent 
with the former results. For the other two pulsars, a reduction in 
the profile width is visible in Fig.Q] The lack of a reliable fidu- 
cial point in the average profiles of these two pulsars prevents 
a robust estimate being made of the retardation and abberation 
time delays, hence estimate of emission heights. 

The drift analysis revealed the large modulation index in all 
three pulsars at this frequency, which are similar to that of the 
Crab puls ar (m — 5 to 8) as found by IWeltevrede et al.l (|2006a) 
and lWeltevrede et al.l(l2007l) at 1380 and 382 MHz, respectively. 
Combined with the values in their work, the modulation index of 
the pulsars considered here increases with decreasing frequency 
indicating the prevalence of bright and spiky e mission. From 
their observations in the 341^4850 MHz range, iKramer et al.l 
d2003l) report giant pulses in PSR Bl 133+16 at 4850 MHz, with 
a majority of them occurring at the phase of the leading com- 
ponent of the average emission profile. They also report individ- 
ual pulses at 1412 MHz that can be distinguished clearly from 
the noise, but have less than 10 times the average pulse energy. 
Furthermore, these pulses were contemporaneous with the giant 
pulse emission at 4850 MHz. If the pulses we detect are related 
to those found at 4850 MHz, then this raises the question of why 
the mechanism responsible for the bright pulses is not active at 
1412 MHz. A possible reason is the unusual spectral behaviour 
of these bright pulses. Alternatively, the average emission may 
have a greater contribution from the normal pulses at 1412 MHz 
than at 4850 MHz or at 116 MHz. Simultaneous observ ations 
spanning even wider frequencies than lKramer et al.l d2~003) could 
help us to answer this question. 

We note that the slopes at the edges of the average emission 
profiles seen in PSRs Bl 112+50 and B0031-07 are indicative 
of the direction of the domin ant drift mode. These p ulsars were 
also identified as drifters bv lWeltevrede et all d2007l) . From their 
work, it can also be seen that several pulsars exhibiting coher- 
ent drifting indicate that one of the leading or trailing edges 
of the average emission profiles is steeper than the other (e.g 
B2012+51, B2310+42, B2016+28 etc.). 

10. Conclusions 

The main conclusion that can be drawn from the pulse energy 
distributions for PSRs Bl 133+16 and B0031-07 is that they 
show a power-law like form in the frequency range observed 
here. PSR Bl 112+50 displays a pulse energy distribution with 
a somewhat steeper power law as a result of several pulses with 
energy >10x(£'). Even though the pulses detected are much nar- 
rower than the average pulse profile, these characteristics alone 
do not qualify these pulsars as giant pulse emitters in the clas- 
sical sense. However, a single pulse analysis for these pulsars 
does reveal large modulation indices, suggesting that they are 
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similar to the "spiky" emiss ion observed in the low frequency 
observations of B0656+14 dWeltevrede et al.ll2006d) . Extending 
this argum ent, this emission may be related to the RRAT-like 
emissions dWeltevrede et a l. 2006b). 

In this work, the spectrum of three pulsars in the 110- 
180 MHz range have been derived. PSR Bl 133+16 has been 
found to exhibit a spectral break in this range, while the spec- 
tra of PSRs Bl 112+50 and B0031-07 display spectra peak 
at frequencies higher than 180 MHz. From our detections of 
narrow and bright pulses, we reported DM changes for PSRs 
Bl 112+50 and Bl 133+16, but no change in the DM of PSR 
B0031 -07. We have detected microstructure in PSRs Bl 133 + 16 
and B0031-07, although with our sensitivity we are unable to 
measure the characteristic intensity variation timescales. In a 
drift analysis of the single pulses in PSR B0031-07, two of 
the three drift modes were detected. The modulation indices 
computed from the pulse stacks of all three pulsars have values 
larger than those at higher frequencies. This shows the presence 
of strong pulse to pulse intensity variations and probably arises 
from the combination of drifting subpulses, single bright pulses 
and emission nulls. While it remains unclear whether these pul- 
sars emit giant pulses, the emission does change considerably as 
one moves to lower radio frequencies. 
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